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TECHNICAL MEMORANDUM X- 64794 

COLUMN DENSITIES RESULTING FROM 
SHUTTLE SUBLIMATOR/EVAPORATOR OPERATION 

INTRODUCTION 


A proposal has been made to dispose of the 100 lb of water (H 2 0)/day 
generated from Shuttle fuel cell operation by a more or less continuous over- 
board dump through a supersonic nozzle after the H 2 0 has been converted into 
gaseous form via the sublimator/ evaporator. Astronomers desire a column 
density of no more than 10 12 molecules/cm 2 so that the induced atmosphere 
will be optically thin at all wavelengths. The purpose of this study is to deter- 
mine the column density resulting from the H a O vent. 

The plume origin is taken to be at (x q , y ) , with the plume axis along 

the x-axis (x q = 13. 7 m , y^ = 4. 6 m), as shown in Figure 1. The line of 

sight S originates at the origin and makes angle y with the x-axis. To sim- 
plify the geometry, only the worst case situation in which the line of sight S 
is in the x-y plane will be considered. 


FIELD POINT 



t 



The plume distribution has been determined semi- empiric ally to have 
the form 


n (cos e) n 

0(r , 6) = C r 2 , (0 < tt /2) 


where n = 6.272. <p{ r, 0) is the molecular current along a radius r mak- 
ing angle 0 with the plume axis. The normalization constant N is deter- 
mined by requiring C 


N„ 


/ <p{r , 0) dA 


hemisphere 


1 


where N T is the total molecular flow rate. 

Since dA ^ = 2irr 2 sin 0d0 , the integral becomes 


N 


2ttN 


T n + 1 ’ 


from which 


N 


j&Liil „ 

2ir T 


ANALYSIS 


The number density at any (r , 0) is given by dividing the current 
<p(r , 0) by the mean radial velocity of the molecules. This velocity is given 
by 


2 



V 


1 003 m/sec 


= / 


2yRT 
(y- l) M 


for H 2 0 at T = 273° K. 
Therefore, 


(n + l) N (cos 0) n 

N(r • e) * < i/2 > 


The column density is found by integrating N(r , 0) along the line of sight 
(LOS) , 


n c = / N ( r * 0 ) ds 
s i 


where Sj is the point at which the LOS S intersects the nozzle plane. Let 
u = s - Sj . From the law of cosines, 

r 2 = r I + u 2 + 2r 1 u sin y 


From the law of sines, 


cos 6 


u 

— cos y 
r 


Finally, 


n = 
c 


(n + 1) N t 

2ffV 


CO 


/ 


(u cos y) n 

n+2 

r 


du 
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Let x = u/r t ; 


n 

c 


(n + 1) N ® 

1 f 

2TTVT! J 


(x cos y) n dx 

(x 2 + 2x sin y + l) 


n+2 

2 


where 




x tan y - v 
o J o 


The integration is easily carried out by the Gauss-Leguerre method 
using eight points. The results are shown in Table 1. From the plot of n 

c 

as a function of angle in Figure 2, it can be seen that the 10 12 molecules/cm 2 is 
violated for look angles smaller than 48 deg from the tail. 


TABLE 1. COLUMN DENSITIES 


r 

(deg) 

*1 

(cm) 

I 

(6) 

n 

c 

(no. /cm 2 ) 

log n 

c 

0 

460 

0.443 

1.95 (14) 

14.29 

10 

701 

0.284 

8.22 (13) 

13.91 

20 

958 

0. 160 

3.39 (13) 

13.53 

30 

1 250 

0.764 

1. 24 (13) 

13.09 

40 

1 609 

0. 0294 

3.70 (12) 

12. 57 

50 

2 092 

0.00850 

8.24 (11) 

11.91 

60 

2 832 

0.00158 

1. 13 (11) 

11. 05 

70 

4 224 

0. 000135 

6.48 (9) 

9.81 

80 

8 229 

1.84 (-6) 

4.54 (7) 

7.65 

85 

16 119 

2.41 (-8) 

3.03 (5) 

5.48 
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Figure 2. Column density versus look angle. 

ORBITAL ACCUMULATION 


Water molecules expelled from the Shuttle will have approximately 8 
km/sec relative velocity to the atmospheric molecules. The collision rate 
will be 


5 





6 = Ntrv 

a a a 


where 


N = 1. 5 x 10 14 molecules/m 3 at 420 km 

a 

cr = 1.04 x 1(T 19 m 2 

a 

v = 8 X 10 3 m/sec 

a 

Q '= 0. 125/sec 


or the time between collisions is 8 sec. 

The mean distance from the spacecraft before collision occurs is 
v/Q « 8 km. Since the water molecules have comparable mass to the atmos- 
pheric molecules, the velocity change from collision will be very large. After 
a few collisions, the expelled molecules will be completely randomized and 
their dispersion will be controlled by a diffusion process. The accumulation 
of H 2 0 along the orbital path will be investigated. 

Consider the H 2 0 trail left by the Shuttle to be initially a line of mole- 
cules along the orbital path with a linear density A molecules/cm. Next, 
consider a cross section perpendicular to this line with thickness df . The 
number of molecules in this volume is Adf. Let cr(r , t) be the area density 
in this cross section (molecules/cm 2 ). The volume density N(r , t) is related 
to (r(r , t) by 

<r(r , t) = N(r , t) di 

Assuming azimuthal symmetry, the molecules will spread in time 
according to the diffusion equation 


9<r (r , t) _ _ , . . 

rr = v . a Vo(r , t) 

ot 
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If a is not a function of r , this equation reduces to 


9£ (r , t) a_ _9_ / 0£ (r , t) 
3t r 9r V 9r 


where a is the diffusion coefficient given by 


a 


A f KT 

mi? 


1 / 7TKT 

,/TN or V 8m 
V a a 


where A f is the mean free path after randomization, K is Boltzmann's con- 
stant, m is the moiecular mass, and T is the temperature. 


The solution to the diffusion equation in polar coordinates has the 

form 


o (r , t) 


C -r 2 /4o!t 
2a t 6 


where C is a normalization constant. This constant is determined by 
requiring 


00 

/ tr (r , t) 27rrdr Adi 
0 


where Adi is the total number of molecules in the section under consideration. 
Carrying out the integration results in 


2ttC = Adi 
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Therefore, 


c (r , t) 


Adi -r 2 /4o!t 
~ — - e 
47rat 


and 


N(r , t) 


A -r 2 /4at 
4wat 


The column density, looking perpendicular to the flight path, is given 
by 


n 

c 


00 

f N(r , t) dr 
0 


A 

47TQit 


/ 

r 

o 


-r 2 /4at 

fcj 


dr 


n 

c 


A 

111 
4ir ia 2 t 2 


The linear density A is just the expulsion rate of molecules N 

T 

divided by the velocity of the spacecraft relative to the atmosphere, or 


A 


1.76 x IQ 22 
8 x 10 5 


2.2 x 10 16 molecules/cm 


Assuming the nozzle is directed such that the molecules are placed behind the 
spacecraft, they will not be reencountered until one orbit later, or t 6 000 
sec. 
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The diffusion coefficient is 


/Z 8.32 x IQ 7 ergs/ mol/deg 270 

_ v8 18 Km/mole )_ 

V2~ (l. 5 x io 8 mol/cm 3 ) (1.04 x 10" 15 cm 2 ) 

= 1.00 x 10 11 cm 2 /sec 


After one period, the column density is 


n 

Ci 


2.2 x 10 


16 


4n/tT (1. 00 X 10 11 6 000) 


1.26 x 18 8 mol/cm 2 


This is the contribution from the first orbit. One orbit later, this will have 
been reduced to 


n 



and the total column density is 



This series does not converge as N — *■ «j ; however, it increases 
slowly. Other factors ultimately place an upper bound on the increase in 
column density. First, the diffusion coefficient increases with altitude. 
Second, oblateness perturbations and earth rotation prevent the spacecraft 
from traversing the same path through the atmosphere on subsequent orbits. 
For a 7-day mission, the number of orbits is ~100 and the series sum is ~18. 
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For 30 days, N » 432 and the series sum is ~40, This is probably a worst 
case estimate because of the preceding consider ati ms, but the maximum effect 
from orbital accumulation for a 30-day mission is 5 X 10® moiecules/cm 2 . 

The molecular impingement resulting from flying through the molecules 
expelled from the previous pass is given by N(0 , 6 000 sec) v or 

£L 



Av 

a_ 

47TOit 


2. 2 x IQ 16 8 X 10 s 
4 tt(i. 0 X 10 11 ) (6 000) 


= 3. 67 x 10® mole/cm z sec . 


Previous orbits will each contribute l/n times this amount. Summing 1/n over 
432 orbits yields 


Q = 2. 2 x 10® molecules/cm 2 sec 

ret 


at the stagnation point. 


MOLECULAR BACKSCATTER 


If the nozzle is oriented so that the expelled molecules are directed in 
front of the spacecraft, the spacecraft will fly through its own cloud as these 
molecules collide with atmospheric molecules. The major contribution will 
come from those molecules that collide fairly close to the spacecraft before 
they become randomized. Therefore, a single collision model is more appro- 
priate than the diffusion model developed previously. 

For simplicity, let the expelled molecules leave the spacecraft with 
velocity v in a parallel beam along the x-axis. Their mean free path rela- 
tive to the spacecraft is x, = v (N <7 v ) _1 . The probability of a molecule 

X r a a a 

colliding with an atmospheric molecule in the interval dx at distance x is 
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P(x) dx = Po(x) [ 1 - Po (dx)] 


where Po(x) is the survival probability after traversing distance x . This is 
given by Poisson statistics as 


Po 



Therefore, 


P(x) dx 


N cr v 
a a a 

v 

r 


exp (-N cr v x/v ) 
a a a r 


dx 


The probability that a molecule will be scattered into solid angle dcu 
is given by 




The differential scattering cross section in the center of mass system is iso- 
tropic and is given by 



Since the velocity of the molecules relative to the spacecraft is small 
compared to the velocity relative to the atmospheric molecules, the expelled 
molecules will be considered at rest in the lab system. Also, the expelled 
molecules will be considered to have the same mass as the atmospheric mole- 
cules. The differential cross section transforms according to 
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2irsin 0 d 6 


2 7rsin 9* d0* 


where 9* is the scattering angle in the c/m system and 6 is the angle at 
which the stationary molecule is scattered relative to the colliding molecule. 
This angle transforms according to 


9 * = tt - 20 


Therefore, 


1 /cur \ _ 4 cos 9 /dtr \ 

o’ \dS2/ (7 \dSV , 

a 0 a e/m 


cos 0 


The scattered velocity is seen from the Figure 3 diagram. 




,♦ v */"» 


Figure 3. Collision geometry for molecules emanating from 
the spacecraft colliding with atmospheric molecules. 


If the two particles have equal mass, v , - v /2. The velocity in 

c/m a J 

the lab system is found by adding w to the x component. 


v v 

a a 

— - — cos 0* 
2 2 


- sin 0* 

Li 
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from which 


v' 2 = (1 - cos 6*) 

z 


Using Q* = r — 29 , 


V 2 = (1 + cos 26) = v 2 cos 2 9 

2 ' a 


or 


v f = v cos 9 
a 


The number of molecules per unit time scattered from dx into solid 
angle cL> is given by 


N T p w dx ^ (fg) 

a 0 


(see Figure 4) 


The molecular flux is this number divided by the area subtended by do , which 
is (r 2 + x 2 ) do; . The contribution to the number density is this flux divided 
by the velocity perpendicular to the area subtended by do? . Therefore, 


dN r (r) 


N_P(x) dx , 

T v /do \ 

(r 2 + x 2 ) v'o- Vd£2y 
3. I 


Inserting the expressions for 



, P(x) dx , and 


v f , gives 
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Figure 4. Contribution to column density from molecules 
backscattered by atmospheric collisions. 


dN»(r) = 


N or v dx cos 0 exp (-N a v x/v ) 
a a a T v a a a r 7 

ir v (r 2 + x 2 ) v cos 6 
r ' a 


The column density contribution from dx is given by 


dn = f dN' (r) dr 
r 

o 


N„,N a dx exp (-N cr v x/v ) 00 

Taa " " " 


a a a r 


r 2 + x 2 


In the iim r -•* 0 , this becomes 
o 
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15 



Therefore, 


^1.5X10 8 (l.04XlO~ 15 ^1.76X10** (11.03) 

n c 2( 1, 00 x io 5 cm/sec) 

n = 1 . 5 x 10 11 mol/cm 2 
c 


The impingement on the spacecraft from the backscatter is given 
approximately by 


Q 


back 


N or N m 
a a T 

4x v R 
r o 


v 

a 


( !• - 108 B)(- 




1.04 




4ir(10 5 cm/sec) (1370 cm) 


= 1.3 x 10 12 mol/cm 2 /sec 


This is two orders of magnitude Less than the impingement of atmospheric 
molecules. 


CONCLUSIONS AND RECOMMENDATIONS 


The expulsion of 100 Ib/day of fuel cell water through a supersonic nozzle 
directed along the x-axis will result in column densities greater than 10 12 
molecules/cm 2 only when viewing at angles less than 48 deg from the tail. 
Smaller viewing angles could be achieved by canting the nozzle slightly down- 
ward. The effect of accumulation of H 2 0 in the orbital path is negligible. Even 
if the H 2 0 plume were directed into the velocity vector, the column density from 
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the backscatter would still be almost an order of magnitude below the 10 12 
molecule/cm 2 criteria, which, incidentally, is two orders of magnitude smaller 
than the column density of the residual atmosphere at 420 km. 

If the spacecraft is randomly oriented, the return flux from atmospheric 
backscatter is ~ 10 12 molecules/cm 2 sec at the stagnation point. This is con- 
siderably less than the flux of atmospheric molecules. If the molecules hit and 
stick to a cryogenic surface, a monolayer (10 15 molecules/cm 2 ) will form in 
10 3 sec. After 7 days, 600 monolayers would deposit. This is 18 micrograms/ 
cm 2 , or 1 800 k . This could be avoided by always orienting the nozzle away 
from the velocity vector or by shielding the cryogenic surface from the velocity 
vector. 
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